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1. The aim of the present work is to obtain generating functions of Chebyshev polynomials 
in three variables. Chebyshev polynomials in several variables associated with the root 
systems of the simple Lie algebras studied intensively in the last few decades n - 0 
Applications can be found for them in different areas of mathematics, as well as in the 
physical investigations. The examples of such applications with the references therein may 
be found in works m-ra. 

Multivariate Chebyshev polynomials are natural generalization of the classical ones. The 
classical Chebyshev polynomials of the first kind Tn{x) are defined by the following formula 
(see for example El) 


Tn{x) = Tn{cos(j)) = COS nf, n = 0,l,.., (1) 

and they satisfy the following three-term recurrence relation 


Tn+l{x) = 2xTn{x) - Tn-l{x), U > 1. 

The classical Chebyshev polynomials of the second kind Un{x) are defined by the formula 


sin [n + 
sin 0 

which is different from o. but they satisfy the same recurrence relation ([2]). The initial 
conditions for the above polynomials 


Unix) 


1)0 


n > 0, X = COS0, 


( 2 ) 


To(x) = 1, Ti(x) = X, Uoix) = 1, Uiix) = 2x (3) 

together with the recurrence relation ([2]) determine them uniquely without reference to (fTj). 

It is known that the function cos nf can be treated as an invariant mean of the exponential 
function on the Weyl group of the Ai algebra root system 

cosn0 = (4) 

Transition to a function which is invariant mean on a Weyl group W of other root system 
leads us to generalization of the classical Chebyshev polynomials to polynomials in several 
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variables Hl-Bl- Let us remind shortly the way in which such hh-invariant functions can be 
constructed. 

Let R{L) be a reduced root system of any simple Lie algebra L. Such a system is a set 
of vectors in d-dimensional Euclidean space supplied with a scalar product The 

system R is uniqually determined by the base of simple roots ctj, i = and a hnite 

group W{R), which is generated by the base reflections. This group is called a Weyl group. 

Generating elements Wi, i = 1, ..,d of the Weyl group act on the simple roots according 
to = —(Uj. The system R is closed under the action of a Weyl group. The elements w 
of W{R) acts on any vector x G by the rule 


2{x,ai) 

WiX = X — - - 

{ai,ai) 

To each root a & R corresponds a coroot dehned by the formula 


(5) 


a 


V 


2a 

(a, a) 


The base of fundamental weights Aj, i = 1,.., d is dual to the coroot base aA, 


1,.., d 


(A.,«J) = 5p (6) 

(we identify the dual space E^* with E^). It is worth noting that the both bases are not 
orthogonal. 

Multivariate Chebyshev polynomials of the hrst kind can be dehned by the following 
function 


(7) 

where the vector n G is given in the fundamental weights base and the vector (p is given 
in the coroot base {a^^} 

d d 

n = ^ riiXi UiE Z, </> = ^ pi G [0,1), (A^, aj) = 6ij. (8) 

^=1 i=l 

Evidently the function $ 11 ( 0 ) is invariant with respect to action of the corresponding Weyl 
group, 


= $n(0), W EW. 

Many of the useful properties of *hn(0) were discussed in the review |5], where these functions 
were called orbit functions. 

In the case of the Lie algebra Ai there is the only simple root a = a'^ (we use the 
standard normalization (a, a) = 2) and there is the only fundamental weight A = a/2. The 
Weyl group W{Ai) consist of the identity e and the idempotent w : wa = —a. Thus, the 
orbit function (jT]) has the form 
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_ ^2ni(n,4)) ^2ni(wn.,4>) _ g7rin</) 

which is identical with (jl]) up to the factor 1/2 and redehnition of phase </>. 

Multivariate Chebyshev polynomials of the second kind can be dehned by the Weyl 
character formula 


^n(0) 


^ detw 

^ det w 
iogW 


<^as 

n+p 

<T)as ’ 
P 


(9) 


where dettc = (—and i{w) is the minimal number of Weyl group generating elements 
Wi required for expressing tc as a product of Wi. p is the Weyl vector which is equal to 
the sum of the fundamental weights. It is essential for our purpose that the numerator and 
denominator in ([2]) can be represented in the following form 


g27ri(i(;m,0) _ E 

uisW,detui=l uicW, det ui=—1 

where m is any vector written in the fundamental weight base. 

2 . Generating functions are a powerful tool as in the theory of classical orthogonal 
polynomials, as in various applications. It is evidently that generating functions are important 
for the studying in multivariate Chebyshev polynomials too. In the work [19] the generating 
functions of Chebyshev polynomials of both kinds associated with the Lie algebra A 2 were 
obtained in explicit form. Recently there was proposed a method of obtaining generating 
function of multivariate Chebyshev polynomials of both kinds [20]. This method uses only 
an explicit form of W - invariant functions ([7j) and is independent from knowledge of recurrent 
relations. In [20] the generating functions of Chebyshev polynomials in two variables associated 
with the Lie algebras C 2 and G 2 were obtained by this method. The method can be summarized 
as follows. 

Since the components of n are integer, scalar product in the function <!)„ (JT]) can be 
represented in the form [wn, (p) = 4>)nk-, and the function itself can be written as 


j'g27ri(wAk,0)^"'fe 

loeW k 



In this formula Mk are the following diagonal matrices 


Mk 


diagle^"^®^'^*'’'^^ g27ri(w2Ak,</)) 


g27ri(w|w|Ak,</))^ 


( 11 ) 


Wi are elements of the Weyl group, and |IT| is the number of elements in the Weyl group W. 

Dehne the matrices Rk = {I\w\ — PkMk)~^, where I\w\ is the identity matrix of size |kL| 
and pk are real parameter. In these notations $n has the form 


* 1 * 11 ( 0 ) ‘l’ni,..,nd(0) 


ni'-.-nJ d“ipi..d°dpjj 


(tr(Rpi--Rpd)) 


Pi=--=Pd =0 
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Let us introduce the new independent variables defined by 


i—1 d—i 

Xi = d>ei(0), ei = (o7^, l,o7^). 

The simple structure of the matrices allows us to express the coefficients *hni,..,nd(0) of 
the function 


ni..nd>0 

in the terms of x*. It is evidently, that the function can be considered as a generating 

function of the multivariate Chebyshev polynomials of the first kind. 

For calculation of generating functions of multivariate Chebyshev polynomials of the 
second kind is required only a minor modification of the above method. The modification 
consists in representation of the functions with alternating signs from the Weyl character 
formula ([H]) as a difference of two terms ([TU]) . This representation allows us to use the 
calculation scheme presented in this section. 

3. The Chebyshev polynomials in three variables were studied in the work [21] out of the 
Lie algebraical context. In particular, in this work were obtained the following recurrence 
relations 


Tk,m+l,n SrTfc 
Tk,m,n+1 

and were calculated some first polynomials using these relations and the natural initial 
conditions To,o,o = 1 , = z, To,i,o = r, To,o,i = 

In what follows we use the representation of the roots of the Lie algebra ^3 by the vectors 
of the Euclidean space with the standard scalar product (.,.)e (see for example [I 8 ]). 
In the formulas below the index E will be dropped. Let e^, i = 1,2, 3,4 be a system of 
orthogonal normalized vectors (cj, Cj) = 6ij of E^. In this case the triple of the fundamental 
roots of is defined as ai = ei — Cj+i, i = 1,2,3. The positive roots consist the set 
Oiij = Ci — Cj, 1 < i < j < 4. All the roots have the norm |ai| = \/2. 

In this representation the fundamental weights are defined by the formulas 

, 3 11 3 111 

+ 2^2 + - ^ei - -e^ - -63 - - 64 , 

^ 1 11111 
A 2 — -«i + 02 + -«3 — 2^1 + 2 ^^ ~ ~ 2 ^^’ 

^ 113 1113 

A3 - + 2^2 + - -Cl + -62 + -63 - -64. 

Weyl vector has the form 


'^k—l,m,n+l 1 l,n5 

T i rri rri rri rri 

fc+ljTn—1,71+1 l,m,n+l -^/c+l,m,n—1 /c—l,m+l,n—1 l,n‘ 

^/c,m+l,n—1 ^/c+l,m—l,n l,m,n* 
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( 13 ) 


p — 2(30^1 + 4^2 + Saa) — Ai + A2 + A3 — + ^2 — 63 — 364). 

The Weyl group W (As) is identified with the group of permutations of the vectors e*, i = 
1,4, with 4! = 24 elements. The generators of this group are the elements tci, W 2 i W 3 which 
permutate the vectors ei and 62 , 62 and 63 , 63 and 64 respectively. 

4. To begin with, let us calculate the orbit function (0 associated with the Lie 

algebra A 3 . Following the scheme presented above, calculate the scalar product {wXk, (f)) 
for all the elements w G fF(A 3 ), taking into account the notation (p = and the 

representation (1T2|) of the roots and fundamental weights by the orthogonal vectors e^. The 
result can be written in the form of diagonal matrices 

Ai = diag(0i, —01 + 02 , 01 , 01 , —01 + 02 , —02 + 03, —01 + 02 , 01 , 01 , —02 + 03, —01 + 02 , 

— 03, 01, —01 + 02, —02 + 03, —02 + 03, —03, —03, —01 + 02, —02 + 03, —03, —03, 

— 03 , —02 + 03 ) 

A 2 = diag( 02 , 02 , 01 — 02 + 03, 02 , —01 + 03, 01 — 02 + 03, 02 , 01 — 02 + 03, 01 — 03, 

— 01 + 03, 01 — 03, —01 + 03, 01 — 02 + 03, —01 + 02 — 03, 01 — 03, —01 + 03, 

— 01 + 02 — 03, 01 — 03, —01 + 02 — 03, —02, —01 + 02 — 03, —02, —02, —02), 

A 3 = diag(03, 03, 03, 02 — 03, 03, 03, 02 — 03, 01 — 02 , 02 — 03, 03, 01 — 02 , — 01 , 01 — 02 , 

02 — 03, 02 — 03, —01, 02 — 03, 01 — 02, —01, —01, —01, 01 — 02, —01, 01 — 02)- 

Using these matrices one can express the function ^n{(p) in the form (HA]) 

$n(0>) = tr(e2"AiAi+n2A2+n3A3))_ (^ 4 ) 

Denoting by the exponential let us introduce the matrices Rk = k = 

1, 2, 3, where pk are real parameters, I 24 is the identity 24 x 24 matrix. In the terms of 
the orbit function has the form 

1 (^(ni+n 2 +n 3 ) 

4>ni,n2,n3(0>) = {tT{Rp,RpM) 

From (IT^ it follows that the function = tii^Rp^Rp^Rp^) expressed in terms of new 

variables is the generating function of three variables Chebyshev polynomials of the first 
kind. Calculation of the function (IT^ for the indices (1, 0, 0), (0,1, 0), (0, 0,1) gives us 

4 * 1 , 0,0 = 62 ;, $ 0 , 1,0 = 4r, $ 0 , 0,1 = 6 A, 

where the new variables are defined by 

^ _ g27ri(/,i g27ri(-</>i+</,2) _j_ g27ri(-(/>2+(/>3) g-27ri03 

P — g 2 m </<2 _|_ g27ri((/>i-(/>2+</i3) _|_ g2vri(-(/,i+(/>3) g27ri(</>i-(/>3) g27ri(-(/>i+(/>2-</i3) _|_ g-27ri02 

^ _ g2m</<3 _|_ g27ri(-</i3+</,2) _j_ g27ri(-(/>2+</,i) 
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In the terms of z,z,r generating fnnction is presented in the following form 


where 


F. 


i^k=3;j=5 

J 2,J,/c = 0 


Pl,P2,P3 


7 7 7 

^Pl ^P2 ^P3 


( 17 ) 


Zp^ = l- zpi + rp\ - zp\ + p{, 

Zp^ = l- rp 2 + {zz - 1 )pI - (z^ + E - 2r)p^ + (zz - l)p^ - rpl + p®, (18) 

Zp3 = 1 - zp3 + rpl - zpl + pI 

It is not difficult to calculate the coefficients Kij^k{z,r,z) from (1T7|) . In the Appendix A 
only those coefficients Ki^j^k are listed whose subscripts satisfy the inequality i < k. This is 
because of the symmetry condition Kij^k{z,r, z) = Kk,j,i{z,r, z). 

Calculation of the orbit function $ni,n 2 ,n 3 with the subscripts (1, 0, 0), (0,1, 0) using the 
generating function flT7|) gives us $ 1 , 0,0 = 62 :, $ 0 , 1,0 = 4r. The reason for appearance of the 
integer factors in these formulas is the presence of stabilizers of weight lattice vectors. By 
definition, stabilizer of n is a subgroup of a Weyl group W consisting of all elements w G W, 
such that wn = n. In the work j5] were obtained the stabilizers in the Weyl group W{A^) 
of all vectors n from the weight lattice. In accordance with this result it is convenient the 
following normalization of the three variables Chebyshev polynomials T„i,n 2 ,n 3 of the first 
kind 


^0,0,0 — 


$ 


0,0,0 


$ 


T ) 77 - 1 , 0,0 rri rri 

ni,0,0 7 5 -^0,772,0 ~A 5 -^0,0,723 


$ 


0,n2,0 


$ 


6 


T 


$ 


771,722,0 


721,72 — 2,0 


Tr 


$ 


0,712,^3 


0,n2,n3 


T 

-L 't 


0,0,ns 


ni,0,713 


711,0,713 


T = $ 

- 7711 , 712,713 ^711,712,713- 


In this normalization all polynomial coefficients are integer. In the Appendix B are listed 
first three variables Chebyshev polynomials of the first kind associated with the Lie algebra 
A 3 . 

5. In this section we calculate the generating function of three variables Chebyshev 
polynomials of the second kind. The function ([9]) can be represented in the form 


a 


<f,as 

^711 + 1,712 + 1,713 + 1 


711,712,713 


^1,1,1 


because of the Weyl vector is the sum of the fundamental weights p = Ai + A 2 + A 3 flT^ . 
Let us rewrite the function $“® as the difference $“■*■ — $““ according to sign of det w and 
introduce the following diagonal matrices 
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^ 1 + = diag ( 0 i , —01 + 02 , —02 + 03 , —01 + 02 , 01 , 01 , —02 + 03 , — 03 , — 03 , —01 + 02 , 

— 03 , —02 + 03 ), 

^2+ = diag(02, —01 + 03, 01 — 02 + 03, 02, 01 — 02 + 03, 01 — 03, —01 + 03, —01 + 02 — 03, 
01 — 03, —01 + 02 — 03, —02, —02), 


^3+ = diag(03, 03, 03, 02 — 03, 01 — 02, 02 — 03, —01, 02 — 03, 01 — 02, —01, —01, 01 — 02), 
^1- = diag( —01 + 02, 01, 01, —02 + 03, 01, —01 + 02, —02 + 03, —01 + 02, —03, —02 + 03, 

- 03, -03), 

^2- = diag(02, 01 — 02 + 03, 02, —01 + 03, 01 — 03, —01 + 03, 01 — 02 + 03, —01 + 02 — 03, 
01 — 03, —02, —01 + 02 — 03, —02), 

^3- = diag(03, 03, 02 — 03, 03, 01 — 02, —01, 01 — 02, 02 — 03, 02 — 03, —01, —01, 01 — 02)- 
In the terms of these matrices the function from nominator of ([9]) has the form 


where Mk± = ea;p(27riAk±). Define the following matrices Rk± = {Ii 2 —PkMk±)~^, k = 1,2, 3, 
where pk are real parameters and /12 is the identity 12 x 12 matrix. In the terms of Mk± 
and Rk± matrices the generating function of three variables Chebyshev polynomials of the 
second kind has the form 


P1,P2,P3 


tr(i?i+i? 2 +-R 3 + — R1-R2-R3-] 


tr(Afi_i_ Af2^ Af3_|_ — Afi_ Af2_ Af3_) 

Using the variables x, r, z introduced previously flT^ we finally obtain the following function 


Pl,P2,P3 

where 

Li{pi,P2,P3) ^ 
L2{PuP2,P3) ^ 
L3{PuP2,P3) ^ 
Li{puP2,P3) ^ 


Li{Pi,P2,P3) + zL2{Pi,P 2,P3) + rL3{pi,P2,P3) + zLi{pi,p2,P3) 


7 7 7 

^Pl^P2^P3 


1 2 242 23 32 , 222,2 I 2 , 4 

1 - P2 - PlP 3 - P1P2P3 - P1P2 - P2P3 + P1P2P3 + P1P2 + P2P3 + P1P2P3, 
pIpIpi + PipI - P2P3 - pIpIp 3 , 

-P 3 pIpi + P2P3P1, 

PspI + P 3 pIpI - PipIpI - P1P2, 


and Zp. were defined by (flSj) . Some first polinomials calculated via this function are listed 
in Appendix C. 
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Appendix A The list of coefficients Kij^k{z,r, z) for the generating function of three 
variables Chebyshev polynomials of the hrst kind z) = Kkj^i{z,r, z)) 


Kooo = 24, 

-^001 = —18^, 

Koo 2 = 12r, 

-^003 = —Qz, 

Kqio = —20r, 
i^oii = 16rf — 6z, 

Koi2 = —12r^ + ^zz + 8, 

^013 = -Gz + 6rz, 

Ko2o = ^6zz - 16, 

Ko 2 i = 20z + 4rz - 

Ko 22 = Wrzz — 6z^ — 6z^ — 8r, 

^023 = -6zz^ + 4rz + 12^:, 


Ko 3 o = 24r - 12^2 - 12^^ 
if’o 3 i = 102 :^^ — 28rf — 6z + 12z^, 

K332 = —Sz'^r — SSz'^r + 16r^ + 16^^ — 16, 
-^033 = Gz^ — 16rz + 4z‘^z — 4z, 

^040 = 82 :^ - 8 , 

-^041 = 4r;2 + 14^ — Sz‘^z, 

Kqa 2 = —6z^ — Oz^ — 4r + Qzrz, 

Koa 3 = lOz + 4^r — 4z‘^z, 

Ko 5 o = —4r, 

Ko 5 i = - 6z + 4rz, 

^052 = -4r^ + 4zz - 8, 

Ko 53 = 12zr - 6z, 


Kioi = 14zz - 8, K 121 = 12z^z^ - 2z^r - 2z^r - 4r^ - 24zz, 

-^102 = —lOzr + 6 z, K 122 = —4z + 3;^^ + 162 :r — Oz^zr + 7zz‘^ + 2 fr^, 

-ft"i 03 = Gz^ — 4r, Ki 23 = 6 z^z — 8 r — 7rzz + 3z^ — 9z^, 

Kin = 12r + 3z^ + 3z^ - 13zzr, K 131 = -lOz^z - lOz^z + 7z^ + 7z^ - 8 r + 29rzz, 

Km = 10 r^ 2 : — IO 2 : — 2z‘^z — 8rz, Km = 8z^r + 122 : — 3^^ + 7zz‘^r — 15z‘^z — 18r‘^z + 4zr, 

Kn 3 = 4r^ — 62 ;^r + 4zz + 8 , Km = —6z* — 4z‘^z‘^ + 2z‘^r + 182 :^r + 4zz — 4r^, 


i^i 4 i = 82 :^^^ — 112 :^ — 42:^r — 4z‘^r + 4r^ + 8 , 
Ku 2 = —2z + 62 :^ — 62 :^fr + 2zr‘^ + 8 : 2 ^^ + 2zr, 
Km = 4z^z — 5zzr — 13z‘^ — 3z'^ + 12 r, 


Km = —4r + 62 :^ + 6 z^ — 4 r 2 :^, 
Ki 52 = 4r‘^z — 2z — 4z‘^z — 4rz, 
Km = —2z^r + 8zz — 8 , 


K 202 = 3K - 4zz - 8 , 
K 203 = -4:zr + 6z, 


K 212 = — 8 r^ + brzz + 82 ;^ + 82 :^ + 12 r, 
K 213 = 4zr'^ — 8 rz — 6 z, 


7r‘^zz — 2z^z^ — hrz^ — hrz^ — 8 r^ — 8 : 22 ; + 16, 

—4z — 4rz‘^z + 42 :^^ + 122:r + 2r‘^z, 

-Or'^z^ - Or'^z^ - 24r + 42:^ + 4z^ + 2 ^ 2 : + 2z^z + 12r^ + I2rzz, 
3z‘^zr — 10r‘^z — 7zz^ + 42:^r — 3z^ + 4rz + 42:, 

24;2^r^ - 22 :^^^ - 7z^r - 7z‘^r + 42:^ - 8, 

-3rz^z + 52 :^^ + 32 :^ + 2r‘^z - 2z, 



K2b2 = — 4 r^ + 12 r + 7 rzz — ?>z^ — 3 ^;^, 
-f^253 = 2r^2: - 2zz^ - Arz + 2z, 

K 303 = 2^^ — 8, 

K313 = -2rzz + 12 r, 


K323 = 2z‘^z'^ - 8zz - 4 r^ + 8zz, 

-^333 = —2zz^ — 2zz^ + 7rzz + 5z^ + 5^^ — 8r, 
i^343 = 22 :^^^ - 2rz^ - 2rz^ + 4r^ - 10;^^ + 8, 
-^353 = —Irzz + 3z^ + 3z^ — 4r. 


Appendix B The list of first three variables Chebyshev polynomials of the first kind 


{Tij^ki 

[z, 

r,z) 

= Tkj,i(z,r,z)) 


Tooo 

= 

1 , 


Tni = 

Tioo 

= 



To 30 = 

^010 

= 

r, 


T400 = 

T200 

= 


- 2 r, 

T310 = 

Two 

= 

rz - 

- 3 z, 

T301 = 

Tioi 

= 

zz - 

~ 4 , 

T220 = 

T020 

= 

2 - 

2 zz + r^, 

T211 = 

T300 

= 

3 z- 

- 3 rz + z^, 

T202 = 

T210 

= 

- 2 r 

+ z^r - zz + 4 , 

Ti 30 = 

T201 

= 

— 2 zr + zz'^ — z, 

T12I = 

T120 

= 

5 z- 

- 2 zz^ + r^z — zr, 

To 40 = 


-3rzz — 3 r + 3z^ + 3z^ + r"^, 
zz - 4z^r + z"^ + 2r^ - 4, 


>^2 




3rzz + z z + 2 r 
irzz - 2r^ + 2z^ 


2z'^z + T^Z 


2J2 


3z 


77:2 


772 „ 


z^zr 
-,2 -2 


.2 -2 


zf- 
4r^ - 
z^zr - 

J2^- 


2z^r 

3z^ 

=2„ 


8zz + 42 ;^r 


4r^ + 4z^r + 32 :^ 2 ;^ 


4r^zz + r^. 


Appendix C The list of hrst three variables Chebyshev polynomials of the second 
kind {Uij^k{z,r,z) = Uk,j,i{z,r, z)) 


f^oo 

= 1 , 

Ul4l 

f^ioo 

= 

Uo30 

Uqio 

= r, 

U400 

U200 

= -r + z"^, 

U310 

Uno 

= zr — z, 

U304 

U401 

= -1 + 

U22O 

U020 

= —2:^ + r^, 

U2II 

U300 

= f — 2zr + z^, 

U202 

U210 

= —+ rz"^ + l — zz, 

U430 

U201 

= —zr + zz‘^ — z. 

U42I 

U120 

= z — zz^ + zr"^ — zr. 

Uo40 


— —2 2 
rzz — z — z ^ 

—r — 2rzz + z^ + z^ + r^, 

-1 + 2zz + - 3rz^ + 

2zr - 2zr^ + rz^ - zz^ + 

- 2rzz + z^z + r - z^, 

r - + z^ - z^z + r^z^, 

z — zr^ + zrz^ — zz^ + zr — z^, 

—zz + — rz"^ — rz^ + z'^z^, 

-2zrz^ + 2 :^ + 2zz^ + zr^ - z - zr^, 

—rz + 3zz — z z + r zz — 1 — rz , 

2rz^ + 2rz^ - 2r^ + 1 - 2zz + 2 ;^^^ - 3r^zz + r^. 
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